Magnetic field correlations in kinematic two-dimensional magnetohydrodynamic 

turbulence 
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The scaling properties of the second order magnetic structure function D% (r) and the corre- 
sponding magnetic correlation function C^ir) are derived for two-dimensional magnetohydrody- 
namic turbulence in the kinematic regime where the ratio of kinetic energy to magnetic energy is 
much larger than one. In this regime the magnetic flux function ip can be treated as a passive scalar 
advected in a two-dimensional turbulent flow. Its structure function D^\r) and the one for the 



magnetic field £>2 (j") are connected by an exact relation. We calculate D'^ l ''{r) and thus D^'^r) 
within geometric measure theory over a wide range of scales r and magnetic Prandtl numbers Pr m - 
The magnetic field correlations follow a r -4 / 3 -scaling law and show an anticorrelation at the be- 
ginning of the Batchelor regime indicative of the formation of strongly filamented current sheets. 
Differences to the full dynamic regime, where the ratio of kinetic to magnetic energies is smaller 
than in the kinematic case, are discussed. 
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I. INTRODUCTION 

Magnetohydrodynamic (MHD) turbulence appears in 
many problems in nuclear fusion research as well as in 
astrophysics and geophysics (see, e.g. Refs. 1 and 2 for 
reviews). Compared to the already complicated hydro- 
dynamic turbulence there are additional degrees of com- 
plexity due to the presence of a magnetic field, the pos- 
sibility of dynamos and the action of the magnetic field 
back on the flow through the Lorentz force. In addition, 
the distribution of energy among the modes is controlled 
through the interaction of three cascade processes, re- 
lated to the three conserved quantities of the dynami- 
cal equations in the inviscid limit: the total energy, the 
cross helicity, and for two dimensions the mean square 
magnetic potential. There is, however, a regime where 
some of the more complicated processes are subdomi- 
nant or absent and where a rather complete analysis can 
be performed: two-dimensional (2-D) magnetohydrody- 
namics with weak magnetic fields. This regime can be 
relevant to astrophysical applications, such as quasi-two- 
dimensional turbulent processes in thin accretion discs. 

I 

The situation we consider is that of a 2-D turbulent 
flow, maintained by some external force, in which a mag- 
netic field is advected. The linearity of the equation for 
the magnetic induction implies that without an external 
seed field no magnetic field will ever be generated. More- 
over, the absence of dynamo action in two dimensions 
P| implies that without permanent driving the magnetic 
field will die out eventually. The action of the magnetic 
field back on the flow can be neglected if the Lorentz force 
is sufficiently weak, i.e., if the magnetic field is sufficiently 
small. Because of the absence of dynamos in 2-D the am- 
plitude of the external driving controls the magnitude of 
the magnetic field. We thus have a regime of kinematic 
evolution of a magnetic field passively advected by the 
turbulent flow which can be realized in any 2-D situation 
for sufficiently small driving of the magnetic field. 



The dynamics of the magnetic field and thus its corre- 
lation functions are nevertheless nontrivial since a tran- 
sient growth and the development of fractal like struc- 
tures on small scales are possible, as evidenced by numer- 
ical simulations. |p| Though the range of applicability 
of this kinematic approach is limited, it shows some in- 
teresting features which we would like to describe here. 

On a technical level, the restriction to a 2-D situation 
allows to connect the evolution of the vector potential 
to that of a scalar field for which detailed scaling pre- 
dictions exist within geometric measure theory. [|7HE| On 
the assumptions of statistical stationarity, homogeneity, 
and isotropy for both the magnetic field and the scalar 
field, geometric measure theory provides a connection be- 
tween the structure functions of the velocity field and the 
scalar field. In particular, the scaling behavior can be re- 
lated and the dependence on the Prandtl number and on 
another parameter that characterizes the strength of the 
driving and appears in 2-D only can be extracted. Using 
an exact relation between the correlation function of the 
scalar and the magnetic field we can transfer these results 
to the magnetohydrodynamic problem and in particular 
discuss the dependence on the magnetic Prandtl number. 
The latter can range from Pr m ~ 10~ 6 in the solar con- 
vection zone to Pr m ~ 10 in the interstellar medium. 
JrO| Within the kinematic regime we are able to go be- 
yond previous applications of geometric measure theory 
to MHD turbulence. QUI 

The outline of the paper is as follows. In sections II 
and III the theoretical background in 2-D magnetohy- 
drodynamics and 2-D passive scalar advection is summa- 
rized. For details on the analysis of the 2-D passive scalar 
within geometric measure theory we refer to our previ- 
ous work. || In section IV both are combined and the 
correlation and structure functions of magnetic fields in 
the kinematic regime are discussed. In section V we com- 
pare with correlation functions and structure functions 
found in full numerical simulations; this highlights fea- 
tures of the interaction between both fields which have 
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to be included when the magnetic field becomes dynam- 
ically relevant. We conclude with a summary in section 
VI. 

p (d t u + (u • V)u) = 

d t B + (u • V)B = 
Vu = 

where v is the kinematic viscosity, r\ is magnetic diffu- 
sivity, p is the thermal pressure, and f u and are the 
external driving for each field. The quantity p is the con- 
stant mass density and po denotes the permeability in a 
vacuum. We estimate typical amplitudes for the veloc- 
ity field and the magnetic field by the root mean square 
values M r . TO .s. = (u 2 ) 1 ^ 2 and i? r . TO . s . = {B 2 ) 1 / 2 where 
(•) denotes a statistical average over an isotropic and 
homogeneous ensemble. If then L is an external length 
scale, the fluid Reynolds number R = u r , m . s L/v and 
the magnetic Reynolds number R m = u r , m . s .L/r] can be 
formed. Their ratio defines the magnetic Prandtl num- 
ber Pr rn = v jr\ = R m /R. The relative size of magnetic 
field and velocity field is measured by the Alfven-Mach 
number M, the ratio of the typical fluid velocity U r . m . s . 



In two dimensions, both the velocity field and the mag- 
netic field can be represented by scalar fields by introduc- 
ing the vorticity 

uje z = Vxu , (7) 
and the magnetic flux function 

B = VipXe z . (8) 
The MHD equations then take on the form 
d t oj + (u • V)w = M- 2 (B • V)j + R~ l V 2 u + f u , (9) 

+ (u • v)tA = i?" 1 W + u , (io) 

with j = — V 2, the z-component of the current density 
and and / w the scalar driving fields . 

The situation we want to consider is one of weak mag- 
netic fields, i.e. very large M, so that the term with the 
Lorentz force density can be ignored. This can be esti- 
mated to be reasonable || if the Mach number exceeds 
the magnetic Reynolds number, i.e. if M > R m . Then 
the equations for the fluid do not depend on the mag- 
netic field any more and the statistically stationary state 
that develops depends only on the strength of the force 
f u and the length scale If on which it acts. In 2-D tur- 
bulence it is the enstrophy that cascades down to the 



II. MHD BASICS 

The equations for a magnetic field B(x, t) in a con- 
ducting fluid moving with velocity u(x, t) are 

(1) 

(2) 
(3) 

to the Alfven velocity va — B r . m . s ./(pop) 1 ^ 2 , 



VA 

/ 2 \ 1/2 

= ( W(K.m.s\ 

\ -^r.m.s. J 

By the last line this is also the square root of the ratio 
between fluid kinetic energy and magnetic field energy. 

With these definitions, the dimensionless form of the 
first two equations reads 

(5) 
(6) 

viscous scales. So if we assume homogeneity and isotropy 
the appropriate smallest scale is the enstrophy dissipa- 
tion length r/u — v 1 ! 2 'e^ 1 ^' ', defined in terms of the en- 
strophy dissipation rate e w = |Va;| 2 ) . 

The quantities we want to study are the second order 
structure function D 2 (r), 

^ B) (r) = (|B(x + r)-B(x)| 2 ), (11) 

and the correlation function C\ (?*)) 

cf ) (r) = (B(x + r).B(x)), (12) 

of the magnetic field. (r) and (r) are related by 

D { 2 B \r)=2{B 2 )-2C i *\r). (13) 

Continuity of the fields on small scales and decay of cor- 
relations on large scales implies the following limiting be- 
havior for small and large distances r: 

° 2 (r >-\ 2(B 2 ) : r^L, (Mj 

whereas it is the other way round for the correlation func- 
tion, 



-VjH (VxB) x B + /wAu + pf, 

Po 

(B • V)u + 77AB + f B 
V -B = 



d t u+{u- V)u = -Vp + M~ 2 (VxB) x B + i?r 1 Au + f u 
dtB + (u • V)B = (B • V)u + R'^AB + f B . 
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C { 2 B \r) 



2(B 2 ) 




* 
L . 



(15) 



Similar quantities can be denned for the velocity field 
and the magnetic flux function ip. The quantity eb = 
?7(|VB| 2 ) is the magnetic energy dissipation rate. Since 
ip obeys in the kinematic limit the equations of a passive 
scalar, we can use our previous results for the scaling of 
a passive scalar in a turbulent fluid to obtain the scaling 
regimes for ip. They can be transfered to the magnetic 
field case by use of the relations |l3) (see Appendix A for 
the derivation) 



4 B) (r) = 2^-AC(r), 



and consequently with ( |l3| ) 



(16) 



(17) 



The quantity e$ = 77(|V-0| 2 ) = v(B 2 ) is the mean flux 
dissipation rate. 



III. 2-D PASSIVE SCALAR ADVECTION 

In the limit of large M the magnetic flux ip does not 
act on the flow field u, so that ([lO]) describes the pas- 
sive advection of a scalar field in the velocity field that 
follows from the hydrodynamic driving. A connection be- 
tween the scaling behavior of the fluid, as contained in 
the velocity structure function, 



D 2 (r) = (|u(x + r)-u(x) 



and its longitudinal part, 

D\\(r) = ^ [ pD 2 ( P )dp, 
r " Jo 



(18) 



(19) 



and the scaling behavior of the second order scalar field 
structure function, 



^(OHh/Hx + r)-^)! 2 }, 



oc r 1 * 2 , 



(20) 



can be found within geometric measure theory. 

One estimates the fractal dimension 8 g of level sets 
ijjQ = ip{x) of the scalar field graph G = (x, -0(x)), and 
obtains in the absence of intermittency corrections, i.e. 
= 2C^\ the upper bound 

2#><4-<# H . (21) 



Building on Ref. 7 and extending the results of Ref. 8 we 
were able to obtain scale resolved and Prandtl number 



dependent upper bounds for Sg 1 ^ , which turned out to be 

rather sharp. || If we assume that the bounds on S g 

are indeed equalities, (|2l| ) gives an upper bound on 
viz. 



- 2 - 2 



d In 7 



In \ / 1 + a Pr n 



V3Pr 



(22) 



In the following we will evaluate the right hand side and 
in particular its dependence on the scaling of the velocity 
field and take this as an estimate for the scaling expo- 
nent Q . A s mentioned, this relics on the assumption 
that several inequalites are in fact equalities and that in- 
termittency effects are absent. A careful analysis of the 
derivation shows that the presence of intermittency in 
the scalar field modifies several approximations and that 
the cumulative effect is difficult to control. It is fairly 
straightforward, however, to allow for intermittency ef- 
fects in the velocity structure function. 

The dimensionless parameter a — / '{(^J^ipr.m.s) 
measures the ratio of the scalar driving rate /ip 2 m s 

1/3 

to the flow advection rate ej and is specific to two- 
dimensional passive scalar advection. For large values of 
a the scalar cannot be advected sufficiently rapidly to the 
small dissipative scales so that a space filling distribution 
with — develops. The distance r is measured in 
units of ijuj and the longitudinal structure function D» (r) 

of the velocity field is measured in units of e 2 / 3 ?? 2 . The 
explicit dependence on the longitudinal velocity structure 
function D\\ allows to substitute numerical or experimen- 
tal findings for the velocity structure function, which very 
often differ from theoretical expectations, and to obtain 
predictions for the corresponding scalar structure func- 
tions. 

As in the previous work Q we use in particular struc- 
ture functions obtained from Fourier transforms of model 
spectra in A:-space and evaluate ( p2| ) numerically. Recent 
experiments on forced two-dimensional turbulenc e [fli 15 
and a number of direct numerical simulations Jl6|-|l8[ sup- 
port the existence of a Kolmogorov-like scaling for the 
energy spectrum, E{k) ~ fc -5 / 3 in the inverse energy in- 
ertial subrange k < kf and a scaling E(k) ~ k~® with 
[3 > 3 for the direct enstrophy incrtial subrange k > kf. 
Recent experiments on 2-D turbulence indicate that the 
velocity increments in both cascade regimes show a non- 
intermittent Gaussian behavior. Jl4],[l9| We do not discuss 
here the role of coherent, large scale vortices which are 
formed by the inverse energy cascade in two-dimensional 
hydrodynamic turbulence and which can effect the spec- 
tra in the enstr ophy s ubrange. This problem is still a mat- 
ter of debate. |^,|l[ We therefore take for our analysis 
the following model spectrum for the amplitudes (|uk| 2 ) 
of the velocity field in a Fourier representation in a peri- 
odic box of size L = 2ir: 
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<|u k | 2 )' 



,-7/3,-11/3,3 

J 1 

k- 3 



k 3 exp 



k-k u 



% < k < fc l5 
k\ < k < kf, 
k f <k<k u = 



(23) 



Note the differences in scaling between the amplitudes 
(|uk| 2 ) and the energy spectrum E(k) due to phase space 
factor, i.e. E{k) ~ k~P~ x corresponds to (lu^l 2 ) ~ k~P. 
For a dependence of the results on (3 the reader is refered 
to our previous work, where we showed that a variation of 
/3 in a range between 2 and 4 left the scaling of Di (r) in 
the enstrophy intertial subrange nearly unchanged. We 
therefore take the exponent (3 — 3 for (|uk| 2 ) [see Eq. 
(p3|)l in all calculations here. 

Assuming stationarity, homogeneity, and isotropy the 
relation between velocity spectrum scaling and the ve- 
locity structure function I?2(r) is given by the volume 
average 

£> 2 (r) = l J |u(x + r)-u(x)| 2 oy, 

= T7 I I u k exp(ik • x)[exp(ik • r) - 1] | 2 dV, 
v Jv k 

= 2^(|u k | 2 )(l-cos(k.r)). (24) 

k 

By averaging over all directions (due to isotropy) in k- 
space the cosine gives rise to the Bessel function Jo(fcr), 



# 2 (r) = 2]T(|u k | 2 )(l-J (fcr))- 



(25) 



The model spectrum ( J23| ) is then substituted and the 
summation in (|25|) is evaluated numerically using a finite 
set of wave numbers, equidistant on a logarithmic scale. 
We checked the independence of our results for several 
wavenumber resolutions . 

The algebraic scaling of D2(r) with respect to r is 
shown in Fig. [j] (thick dash-dotted line). For the dissi- 
pation scales r < r/ w it follows a r 2 -dependence in cor- 
respondence with the Taylor expansion. This quadratic 
scaling with respect to r continues in the enstrophy iner- 
tial subrange for scales iq u <r < If. For distances If < r 
the energy inertial subrange with a r 2 / 3 -scaling sets in. 
Finally, a saturation to a constant value due to finite size 
effects is observed. 



IV. 2-D KINEMATIC MHD TURBULENCE 



Using the energy spectra ([23|) we can now evaluate the 
passive scalar structure and correlation function expo- 
nents and, via ([l6[), the corresponding behavior for the 
magnetic fields. The results for a large range of length 



scales and magnetic Prandtl numbers are shown in Figs, 
[j] and ^|. The structure function is plotted in units of 

is obtained by numerical inte- 



-i/r 



and log 10 D 2 



gration of (r) over log 10 r. We find in both figures 
the following scaling regimes. For small scales r <r\ u} the 
fields are smooth and one finds D^\r) = {e^/2r))r 2 and 



thus £. 



2. The quadratic scaling continues into the 



enstrophy inertial subrange r > r]^. For sufficiently large 
Pr m the third term under the square root in Eq. ( j22] ) 
gives the dominant contributions and the scaling expo- 
nent changes to Q = 0; this is the viscous-convective 
Batchelor regime [E2[ . The flux tp is then advected chaot- 
ically in a velocity field that is still smooth on these 



scales, i.e. Dm 



f 2 and the flux contours are stretched 



and twisted by the fluid motion and form filamented cur- 
rent sheets. The r 2 / 3 -scaling regime for larger r and for 
both the velocity field and the magnetic flux function 
is connected with the inverse cascade of the underlying 
fluid turbulence; this is the inertial-convective Obukhov- 
Corrsin regime p3[ . The final saturation to a constant 
value for very large separations is due to the finite sys- 
tem size L. 

The scaling of the structure function D 2 (r) as well 
as of the correlation function C\ (r) follow from Eqns. 



( |16| ) and (17), respectively. The results for the structure 

function of the magnetic field D% (r) in units of e^^" 1 
are shown in Fig. y (thick solid lines). Clearly, for scales 
r <r\ LL> the structure function shows again the Taylor ex- 
pansion behavior D^\r) — (es/2r7)r 2 as stated in Eq. 
(|l4|). The larger the scales the more dominant the con- 
stant first term in Eq. ([if]) , which is proportional to the 
magnetic Prandtl number, leading to a saturation at a 

( B) 

constant value D 2 (oo) = 2e^/rj for large distances, i.e. 



where ~ r u . The crossover between both scaling 

regimes is relatively sharp and shifted towards smaller 
r for growing values of Pr m . We could not observe an 
intermediate algebraic scaling between both ranges. 
The spatial correlations in the magnetic field are sub- 

( B) 

dominant in D\ (r) for scales 

(B) 



> 



rjuj, but show up 

in the correlation function C^'ir): an algebraic scaling 



D^\r) ~ r a corresponds to C 2 D> {r) ~ r a ~ 2 . The results 
for the correlation function, derived using ([l7]), are shown 
in Fig. H for different values of Pr rn and a. In the loga- 
rithmic plot the algebraic decay of the correlation for all 
discussed parameter sets is clearly visible in particular, 



a r 



" 4 / 3 -decay is found in regions where D^\r) 



r 2/3 



(see also Fig. ||). The correlations decay very fast when 
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the magnetic flux is advected passively with the turbu- 
lent velocity field. 

The results show an anticorrelation of the magnetic 
fields, i.e. C\ < 0. It is observed near the crossover from 
', i.e. near the transition from the 
viscous regime to the Batchelor plateau. This follows also 
analytically within the Batchelor parametrization J24| of 

(r) by setting 



~ r 2 to 



(26) 



which approximates such a crossover very well. Equation 
( |l7| ) then gives 



-d r (rd r D^'(r)) 



1 - a2r 2 
: (1 + a 2 r 2 f 



(27) 



With 0,2 > and r > we find negative values of (27) for 



r > (a^) -1 / 2 . It can be shown that a crossover of £>2 (r) 
from r a to r 13 with < (3 < a < 2 gives negative values 
in (^t]) and thus anticorrelation of the magnetic field at 
all when in addition the constraint P < a/5 is satisfied. 
This is consistent with our results for extremely small 
magnetic Prandtl numbers where no anticorrelation is 
observed. Here, the quadratic scaling with r changes di- 
rectly to an r 2 / 3 -scaling (see Figs. || and [I]), and the 
intermediate Batchelor plateau is not present. 

The anticorrelation in the normalized correlation func- 
tion is magnified in Fig. || for values of Pr m reach- 
ing from 1CP 5 to 10 3 . The negative contributions van- 
ish for sufficiently low Pr m where the Batchelor plateau 
is supressed (see the cases Pr. m = 1CP 4 , 10~ 5 ). On 
the other hand we observe a nearly constant value 

min ^Cj 5 ^ — — 0.0375e^,/?7 for the larger values of the 
magnetic Prandtl number. 

The origin of this anticorrelation presumably lies in 
the advection of ?/> in turbulent flow structures and the 
formation of strongly filamented flux sheets (and thus 
current sheets) in which the magnetic energy is stored. 
The minimal sheet width des can be determined by the 
balance between the advection and the diffusion term of 
©, giving 



±- 

>) 2 
J cs 



(28) 



1 /3 

By taking Ui/l — ej as the typical strain rate we end 
up with 



S cs ~ ^Pr- 1 / 2 



(29) 



This Pr m scaling is also found in the position of the min- 
imum, so that the anticorrelation is due to layers with 
opposite orientations of magnetic fields. The saturation 
of the anticorrelation in Fig. || suggests that the ratio of 
regions with perfect parallel alingment and deviations, 
say due to folds or modulations in thickness, stays the 
same, independent of magnetic Prandtl number. 

The relation ( p9[ ) can also be derived by means of (|2^) . 
Assuming the subdominance of the aPr m r 2 term, we can 
ask when the last term under the square root in ( p2| ) will 
exceed unity. By requiring a smooth flow in the viscous 



subrange, i.e. setting Du 
dependence as in (129). 



we get the same Pr„ 



V. FULL MHD BEYOND THE KINEMATIC 
REGIME 

In the following the results of the kinematic approach 
are compared with direct numerical simulations (2^] as 
well as integrations of spectral transfer equations within 
closure models [^6| , both at Pr m — 1 and with a kinetic 
energy about three orders of magnitude larger than the 
magnetic energy. This provides insights into the feedback 
of the magnetic field on the flow via the Lorentz force 
density. The question is whether the phenomena seen in 
the kinematic case, in particular the anticorrelation, will 
still be present. 

Both numerical experiments support the existence of 
a Kolmogorov-type scaling for the spectrum of the mean 
square magnetic potential, A{k) ~ k~ 7/>3 for fc < kf, in 
the inverse mean square flux cascade range. The Alfven 
effect causing an equipartion between kinetic and mag- 
netic energy is subdominant. The spectral closure sug- 
gests A(k) ~ fc~ 7 / 2 for the direct energy cascade range 
k > kf. We start with a corresponding model spectrum 
for the amplitudes (IV'kl 2 ) (see Fig. [|(a)) 



, -7/6,-13/3,3 

f 1 



-7/6 



k -6/2 



exp 



fc -4/3 
k -6/2 
2 



k — kd 



x < k < fc i> 
ki < k < kf, 
kf < k < k d , 

kd < k , 



(30) 



where kd ~ r\ d is the wave number at the viscous cut 
off. The structure function is calculated as in section III: 

^ ) (r) = 2^(|^ k | 2 )(l-Jo(fcr)). (31) 

k 



The resulting structure function is shown in panel (b) of 
Fig. The inset gives (2 vs - r measured in units of 
rjd- Consequently D^ yT) has to be taken now in units 
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of e^e 1 l 2 v 1 l 2 . The r 2 -scaling for r < rjd is followed by 

a range with (3 ~ | for r > 77^, which seems to be 
connected with the extended inverse cascade of the mean 
square magnetic potential. Again structure function and 
correlation function are calculated by means of @ and 
(|l7|), respectively Consistency was checked by calculat- 
ing (r) directly via (|Bk| 2 ) ~ k 2 (\ipk\ 2 } using a rela- 



( B) 

tion equivalent to (31). The function D2 (r) is shown in 
Fig. |3| as the dash-dotted line. We see that it reaches the 
constant saturation state for scales r larger than in the 
kinematic case. The differences between the kinematic 
and full dynamic situation are more pronounced in the 
correlation functions [see Fig. |(c)]. The correlations de- 
cay slower with respect to r, following now a r~ 2 / 3 -law. 

Besides the different scalings we note that the anticor- 
relation has disappeared. In the full dynamic case the 
magnetic field itself would prohibit the build-up of elon- 
gated current sheets (flux sheets) due to its feedback via 
the Lorentz force density 



fi= V(^)+-(B-V)B. 

\2(j, J [i 



(32) 



The current sheets, if formed, become sensitive to sev- 
eral resistive instabilities such as the spontaneous growth 
of the tearing mode, p7|-p9t where the stored magnetic 
energy can be released into the plasma flow by means 
of magnetic reconnection. In this process the Lorentz 
force density causes the plasma acceleration due to re- 
connected magnetic flux and prevents a further steep- 
ening of flux up to the dissipative scales. Thus in ad- 
dition to the condition M > R m the Reynolds number 
R has to be sufficiently low to avoid resistive MHD in- 
stabilities. It was found that the onset of the tearing 
mode instability is determined by the Hartmann num- 
ber Ha = y/R^R = RVPr^ and the width S C s- pOpl 



VI. SUMMARY 

Our findings for the magnetic structure function as 
well as the corresponding correlation function in a two- 
dimensional MHD system can be summarized as follows. 

(1) For the kinematic regime of two-dimensional MHD 
turbulence (M ^> 1), the magnetic flux function ip can be 
treated as a passive scalar advected in a turbulent flow. 
By means of the geometric scaling theory the second or- 
der structure function ( r ) was calculated over a wide 
range of scales r and of magnetic Prandtl numbers Pr m . 
For this case a rather complete analysis is now possible 
by means of Eq. (p2[). 

(2) The second order magnetic structure function 
Z?2 (r) and the corresponding magnetic correlation 

( B) 

function C2 (r) were calculated by exact analytical rela- 
tions from {?). The larger Pr m the smaller the scale 
r where the structure function reaches the constant satu- 
ration. The correlation function decays fast with r~ 4 / 3 



and shows an anticorrelation. The latter is connected 
with the transition to the Batchelor regime and can be as- 
sociated with an attempt of the system to concentrate its 
magnetic energy in strongly filamcnted structures. This 
process is limited by the finite resistivity rj. The anticor- 
relation minimum was found to be independent of the 
magnetic Prandtl number. 

(3) The results for the kinematic approach were com- 
pared to results in the dynamic regime M < R m . The 
magnetic field correlations decay with r~ 2 / 3 , i.e. slower 
than for the kinematic case. This can be explained by 
the feedback of the magnetic field on the flow which pro- 
hibits the build-up of strongly filamented, elongated cur- 
rent sheets on small scales. Therefore no anticorrelation 
was detected. 

While the discussion in this paper was limited to the 
2-D case, a weak field regime without dynamo can also 
be identified in three dimensions. In the experiments of 
Odier et al. j32| a scaling of the magnetic field was ob- 
served that is very close to what one expects for a pas- 
sive scalar. Thus the stretching of the magnetic field does 
not seem to have a pronounced effect on the structure 
function. Further investigations of this point seem worth- 
while. 
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APPENDIX A: DERIVATION OF EQUATION (16) 

( B) 

The second order structure function D\ (r) was de- 
composed in Eq. (|l^) into (the sum convention is used) 



Di B \r) = 2(B? 



-2Ci B \r), 

(B) 



2i?n0)-2i?nr), 



(Al) 



where R\j (r) = (Bi(x)B 3 (x + r)) is the static magnetic 
field correlation tensor. A detailed analysis of the proper- 

( B) 

ties of R\: ° ; (r) 

can be found in Oughton et al. 



we apply tijhUir, 



{Sji5 k 



Next 



SjmSki) for B = VxA in 



the autocorrelated part and get 

= (dx-Am) 2 — {d x -A m d Xm Aj). 



The statistical average gives 

(S 2 ) = ((d Xj A m ) 2 ) = 



(A 



(A2) 



(A3) 



where the second term vanishes due to homogeneity. The 
trace of the magnetic field correlation tensor can be de- 
composed in a similar way by using pi] 



G 



(S i (x)B i (x + r)) 



-e lj keumd rj d ri R ( k f l l{r), 
-d 2 R {A) (r) 



i 



-ADi A \r). 

For the derivation the Coloumb gauge V ■ A 
be taken. Finally we get in 3-D 



D 



IB) 



, (r) = 2^ 



AD {A \r), 



(A4) 
has to 

(A5) 



and it follows in 2-D with A = (0, 0, tp(x, y)) Eq. |[|) 



D 



(B) 



2 (r)as2 a 



A#(r) 



(A6) 



The dimensionless form of (|l|), D [ 2 B) = 2Pr m -AD^\ 
was solved numerically over a large range of scales r. We 
substitute the Laplacian in the following way 



AD 



MO 



D 



WO 



[(dinrg) 2 + d? nr g] 



(A7) 



r- 

Here (/(lnr) = hiD^ir) was set. The function g is ap- 
proximated with cubic splines in a least square fit using 
the Numerical Algorithms Group (NAG) library. 
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FIG. 1. Velocity structure function D2(r) (dash-dotted 
line) and the scaling exponent y ( r ) °f the struc- 
ture function (r) for three values of the prefactor 

a — £ipeZ 1 ^ 3 Tpr.rn.s. an d three values of the magnetic Prandtl 
number. The solid line stands for a = 10 -4 , the dotted for 
a = 10 -3 , and the dashed line for a = 10 -2 . For small dis- 
tances r the a-dependence is subdominant and thus the curves 
for the same magnetic Prandtl number coincide. 
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FIG. 3. Normalized structure function 
D 2 B \r)/(€^u~ 1 Pr m ) vs. r* = r/rj w for a as in Fig. |l|. The 
curves for the different values of a, but the same Pr m , coin- 
cide. Additionally, D ( 2 B \r)/( e^v 1 ) vs. r* = r jr\d calculated 
from direct numerical simulations (DNS) at Pr m — 1 is plot- 
ted as a dash-dotted line. 
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FIG. 2. Structure function D 2 (r) with respect to r 
for a as in Fig. [I]. Additionally the structure function for 
Pr m = 10~ 5 for a = 10 -4 is shown which has no Batchelor 

plateau. All linestyles are chosen as in Fig. fjl. FIG. 4. Normalized magnetic correlation function 

C^ S) (r)/(e^- 1 Pr m ) vs. r/jfc for Pr m = 10 A . For Pr m = 1 
the correlation functions are shown for the three values of a 
as in Fig. hi using the same linestyles as therein. 
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FIG. 5. Anticorrelated part of Cj (r) where the magnetic 
Prandtl number Pr m decreases by a factor 10 from 10 3 for 
the left most curve to 1CP 5 at the right end. The parameter 
a was set to 1CF 4 . 
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FIG. 6. (a): Mean square potential spectrum (|V'k| 2 ) vs. k 
resulting from full MHD simulation results. The extension of 
the two cascade ranges is identical to those of the velocity in- 
put spectrum for the kinematic calculations, (b): Correspond- 
ing structure function D^\r) vs. r plotted with the same 
linestyle as (a) and compared with its kinematic counterpart 
for a — 10 -4 . The inset shows the corresponding ^\r). (c): 
Both magnetic correlation functions are plotted correspond- 
ing to (b). All quantities are expressed in characteristical units 
and Pr m = 1 was taken. 
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